Prom cluster to solid - the variational cluster approximation applied to NiO 
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The variational cluster approximation is applied to the calculation of the single particle spectral 
function of NiO. Trial self energies and the numerical value of the Luttinger-Ward functional are 
obtained by exact diagonalization of iViOe-clusters and the single particle parameters of the clusters 
serve as variational parameters to obtain a stationary point of the grand potential of the lattice 
system. Good agreement with experiment is obtained. 
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The theoretical description of compounds containing 
partially filled 3d, 4/ or 5/ shells is a much-studied prob- 
lem in solid-state theory. Due to the small spatial ex- 
tent of these shells the Coulomb repulsion between con- 
duction electrons in these compounds becomes unusu- 
ally strong and approximations which rely on a map- 
ping of the physical electron system onto one of fictious 
free particles in a suitably constructed effective poten- 
tial cannot even qualitatively describe the resulting state. 
Starting with the work of Hubbard [l| a variety of the- 
oretical methods have been invented to deal with this 
problem @, & H H IE 0, II H • Ma J or progress towards a 
quantitative description of 3d transition metal (TM) ox- 
ides has been made by the cluster method initiated by Fu- 
jimori and MinamijlollTH. This takes the opposite point 
of view as compared to band theory, namely to abandon 
translational invariance and instead treat exactly - by 
means of atomic multiplet theory [12, HH - the Coulomb 
interaction in the 3e2-shell of a TM-ion in an octahedral 
'cage' of nearest-neighbor oxygen atoms. Recently ideas 
have been put forward to broaden the spectra of finite 
clusters into bands [l4l US [IB]- ln particular, building 
on field-theoretical work of Luttinger and Ward[l7j who 
showed that the grand canonical potential f2 of an in- 
teracting Fermion system is stationary with respect to 
variations of the electronic self-energy S, Potthoff has 
recently proposed[l||] the variational cluster approxima- 
tion (VCA) where trial self energies are generated numer- 
ically by solving finite clusters and used in a variational 
scheme for J7. So far the VCA has been applied mainly 
to simplified systems such as the single-band Hubbard- 
model[H, HH but the success of the cluster method for 
TM-oxides suggests to apply the VCA also to a realis- 
tic model for TM-oxides thereby using the octahedral 
clusters discussed above to generate self-energies. Here 
we outline such a calculation for the frequently studied 
compound NiO. Using clusters containing just a single 
TM-ion implies that the self-energy is site-diagonal, i.e. 
fc-independcnt. The corresponding approximation thus 
is similar to the dynamical mean-field calculations which 
have recently been applied to a variety of compounds [201] . 
We start by defining the Hamiltonian which describes 
the NiO lattice and denote by d\ a a an operator 
which creates a spin-cr electron in the ci-orbital a € 
{xy, xz, yz, . . . } on metal site i and Pj,x,a annihilates an 



electron in p-orbital A S {x, y, z} on oxygen site j. Tak- 
ing the 02p-level energy e p as the zero of energy the 
single-particle terms read 
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The Hamiltonian also contains terms which describe hy- 
bridization between next-nearest neighbors, i.e. 02p- 
02p and Ni3d-Ni3(i. The numerical values of the differ- 
ent hopping integrals have been obtained by an LCAO- 
fit to a paramagnetic LDA band structure of NiO[2lj. 
They are quite similar to those used by Fujimori and 
Minamifioj and van Elp et a/.[ll[. Also included - via 
the e Q - is a CEF-splitting of lODq = 0.7eV^, [Tl|. 
The Coulomb interaction within the d-shell reads 
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where we have suppressed the site label i and A = (a, a). 
The matrix elements V^'^ can be expressed [lU, [l3|] 
in terms of the 3 Racah-parameters, A, B and C . To 
avoid an 'implicit' interaction between electrons these 
are assumed to be independent of the d-shell occupation. 
B and C can be estimated from atomic Hartree-Fock 
wave functions and the values used here, B = 0.13eV 
and C = 0.6eV are standard ones 0, El. The pa- 
rameter A is reduced substantialy by solid state screen- 
ing and is usually adjusted to match experiment - here 
we do the same and choose A — 9eV. The 'Hubbard 
U'= E(d 9 ) + E(d 7 ) - 2E(d 8 ) - calculated from ground 
state energies of free ions without CEF-splitting - isfl3l] 
U = A + B = Q.leV. Previous estimates range from 
U = 6.7eV - by fit of cluster spectra to experime nt fToL Hl| 
- to U = 8eV by density functional calculations [2J]. The 
present value thus is somewhat large which will be dis- 
cussed below in more detail. Next, we choose - again 
by adjusting to experiment - — e p = —68.5 eV so 
that the charge transfer energy A = E(d 9 L) — E(d s ) 
is[I3 A = e d - e p + 8A ~ 6B + 7C = 6.9eV. These val- 
ues are consistent with the notion that NiO is a charge 
transfer insluator but close to the intermediate regime 
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of the Zaanen-Sawatzkv- Allen diagram(23j. Finally, any 
Coulomb interaction between electrons which are not in 
the same 7Vi3d-shell is neglected. 

The VCA is based on an expression for the grand po- 
tential f2 of an interacting many-Fermion system due to 
Luttinger and Ward [17]. In a multi-band system where 
the Green's function G(k, uj), the nonintcracting kinetic 
energy t(k) and the self-energy S(k, uj) for given energy 
uj and momentum k are matrices of dimension 2n x 2n, 
with n the number of orbitals/unit cell, it reads [13] 

il = -iy e ^ 0+ ln deif-G^fk,^)) +F[S] (3) 

where uj v — (21 + (3 are the Matsubara frequencies, 

G _1 (k,a;) =w + //-t(k)-S(k,w) (4) 

and the functional F[S\ is the Legendre transform [l8j 
of the Luttinger- Ward functional $[G]. A definition of 
$[G] in terms of Feynman diagrams was given by Lut- 
tinger and Ward(i~7l], a nonperturbative derivation has 
recently been given by Potthoff [2{| . f2 is stationary with 
respect to variations of the self-energy [l7| 



but a prohibitive obstacle in exploiting |5} in a varia- 
tional scheme for X is the evaluation of F[E] for a given 
'trial E\ Potthoff has suggested 18] to restrict the do- 
main of S to 'cluster representable' ones, i.e. exact self- 
energies of finite clusters, for which F[S\ can be deter- 
mined numerically from The key observation [l8| is 
that $[G] and hence its Legendre transform F[S\ have 
no explicit dependence on the single-particle terms of H 
whence F[S\ is the same functional of X for any two 
systems with the same interaction part of the Hamilto- 
nian. Under the assumption that only interaction lines 
connecting orbitals in the same d-shell are relevant we 
can therefore use a numerically soluble system of discon- 
nected finite clusters - the so-called reference system - to 
generate trial self-energies E(u>) together with their ex- 
act -F[X]. More precisely we choose a reference system 
where each Ni3d orbital d a is coupled to one 'ligand' or- 
bital L a with these ligands in turn decoupled from each 
other and the interaction within the c?-shell given by 
The reference system thus is equivalent to an array of 
non-overlapping identical Ni06 clusters where each lig- 
and L a corresponds to the unique linear combination of 
02p orbitals on the six nearest O neighbors of a given 
Ni atom which hybridizes with the Ni3d a orbital. After 
numerical diagonalization of the cluster Hamiltonian for 
all possible electron numbers we obtain the grand poten- 
tial O and Green's function G(u>) of the cluster whence 
equations (|4]) and ([3]) give and F[E]. Next we in- 

sert the and F[E] so obtained into equations ^ 

and Q for the lattice system and obtain an approximate 
G(k, uj) and fl for the infinite system. Variation of S(oj) 



is performed by varying the single-electron parameters - 
such as hybridization integrals or site-energies - of the 
reference system and the best approximation to S(w) is 
obtained by demanding that Q be stationary with respect 
to such variations. We write the single-particle Hamilto- 
nian for a N1L5 cluster as 

H 8ingle = V ( a ) ( di,, L a ,a + H.c. ) 

a,cr 
a,cr 

and have the following varational parameters: 

1) The hopping integrals V(a). Since the ground state of 
rf 8 in cubic symmetry has the configuration t\ g e^ g , V(t 2g ) 
connects mainly completely occupied orbitals. To sim- 
plify the problem we therefore discard the three ^g-hke 
ligands alltogether and write the remaining V(e g ) = Xt 
where t — \/2>(pda) is the hop pin g integral in the cluster 
calculation for a NiOfi-cluster|10l|. 

2) The site energy of the e 9 -like ligands e(e g ). 

3) The site energies E(a) of the e g and t 2g -\ike d-orbitals. 
By numerical scan a set of these 4 parameters where 
n is stationary can be found for each temperature T. 
For temperatures between 100 — lOOOif the expression 
fi(T) = ^0 — fcsTlog(3) - exptected for a gapped, para- 
magnetic spin-1 system - gives an excellent fit to the cal- 
culated f2(T). The constant fio thereby is the ground 
state expectation value (H — fiN) . Since the VCA also 
gives the lattice Green's function G(k, u>), (H — fiN) can 
alternatively be computed from its th and I s * moments. 
There is no reason why these two results should agree 
- the formalism of the VCA would not necessitate this. 
Still the two values agree quite well — 418.649 eV/unit 
cell versus —418.684 eU/unit cell - indicating that the 
VCA is a quite 'intrisically consistent' approximation. 
Next, we proceed to a comparison of G(k, uj) to exper- 
iment. Figure [T] compares fe-integrated spectral densi- 
ties calculated at room temperature (300 Kelvin) to an- 
gle integrated valence band photoemission spectra taken 
by Oh et al. [26[ at two different photon energies. At 
hv = IbOeV the experimental spectrum essentially re- 
sembles the d-like spectral density, whereas at hv = 67eV 
the states at the valence band top are anti-resonantly 
supressed - whence 02p-derived features become more 
clearly visible - whereas the 'satellite' at — lOeV is reso- 
nantly enhanced. Figure [2] compares the fc-resolved spec- 
tral function for momenta along (100) (r — » X) to the 
experimental band dispersion by Shen et aL[2?|. 

The spectral density has gap of approximately 4eV 
around the chemical potential. This is consistent with 
experiment 28] but has of course been achieved by the 
choice of A and A. At the top of the photoemission 
spectrum E < there is a high-intensity band complex 
at binding energies between w —4 eV and » —2 eV, 
which was shown to consist of several sub-peaks by Shen 
et aL[13]. These authors did not actually resolve the 
dispersion of the individual sub-peaks although the data 
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(a) 02p-like 

Ni3d-like - 
XPS, hv=150eV • 




(b) 02p-like 

Ni3d-like - 
XPS, hv=67 eV • 
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FIG. 1: Single particle spectral densities obtained by VCA 
compared to valence band photoemission data (XPS). 



seem to indicate a weak overall 'upward' dispersion as one 
moves r — > X which would be consistent with theory. It 
also has to be kept in mind that the calculation has been 
performed for the paramagnetic phase whereas the ex- 
periment was done below the Nccl temperature and thus 
in the antiferromagnetically ordered phase. This may 
have an impact on the dispersive features. Proceeding 
to more negative binding energy both the experimental 
band structure and the theoretical spectra show a gap of 
~ leV and then a group of dispersionless bands between 
—4.5eV and — 6.5eV. In the angle-integrated spectrum, 
FigureQJi, the topmost of these bands produces the shoul- 
der at — 4.5eV. The spectrum in Figure []J> shows a peak 
at —6eV which originates from states with essentially 
pure 02p character [26j . The corresponding peak in the 
theoretical spectrum originates from a saddle point sin- 
gularity of the upper 02p derived band at the X-point. 
Comparison of the angle integrated spectrum, the exper- 
imental dispersion and an LDA band structure makes 
this a plausible explanation. The sole strongly dispersive 
feature in the spectrum, namely an (92p-derived band at 
binding energies between —6eV — > —9eV is again well 
reproduced by theory. Finally the 'satellite' at bind- 
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FIG. 2: Top: Dispersion of experimental bands measured 
by ARPES 27]. Bottom: fe-dependent spectral function for 
momenta along V ~ X. Lorentzian broadening 0.1 eV, d-like 
weight is multiplied by factor of 5. 



ing energies —9eV — > — IheV consists of at least two 
sub-peaks as can be seen in Figure |T|b and also in the 
ARPES data. Theory predicts several sub-peaks in the 
satellite but these may not have been resolved in exper- 
iment due to the strong broadening of the satellite. By 
and large we may say that there is essentially a one-to- 
one correspondence between theory and the measured 
bands (the k-integrated spectra also agree roughly with 
a recent LDA+DMFT(QMC)-calculation@). Finally we 
mention the values of the ri-shell occupation rid = 8.16 
and the expectation value of the d — p hybridization 
(H pd ) = -2.30 eV/unit cell (at T = 300 K with neg- 
ligible T-dependence). To put these in perspective we 
note that perturbation theory for a NiOg cluster with un- 
renormalized parameters gives rid = 8 + 2 (t/A) 2 = 8.21 
and (H pd ) = -4P/A = -2.93eV. 
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An issue that requires more detailed discussion is the 
spectral weight of the satellite: the fraction of d-like 
weight below the'gap' around — 7eV is w — 61%. A 
recent LDA+DMFT calculation Q gave w = 45% with 
U = 8eV. Van Elp et a/.[ll| estimated w = 30% with 
U = 6.7eV. As shown in Ref.flil] the spectral weight of 
the satellite is directly related to the choice of U - larger 
U produces a more intense satellite. The relatively large 
value of U = 9.1eV in the present work was necessary to 
get the satellite below the bottom of the 02p-band - as 
suggested by ARPES[13]. This raises the question why 
we need a larger value of U as compared to the 6.7eV 
in the cluster calculations [13, HH • ^ turns out that the 
reason is the downward renormalization of the hopping 
integral V(e g ) in (J6j> . The value of A at the stationary 
point is 0.812 (at 300K with negligible T-dependence) 
and even such a small reduction leads to an appreciable 
upward shift of the satellite due to reduced level repulsion 
between satellite and valence band top (see e.g. Figure 11 
of Ref.[ll|) - which must be compensated by a larger U. 
On the other hand there is a simple physical argument for 
the value of A < 1: in a NiOg cluster, the mixing strength 
between (say) a Ni3d x 2_ y 2 orbital and the bonding com- 
bination of 02p x orbitals on the two nearest neighbors in 
x-direction is t/y/2. In a NiO lattice the matrix element 
between the Bloch states of momentum k formed from 
the Ni3d x 2_ y 2 orbital and the 02p2;-orbital is itshi(k x a) 
with a the Ni-0 distance. The Brillouin zone average of 
|sin(fca;a)| is 0.63 = 0.89/-\/2 so that the value of A at 
least partly reflects the smaller average d — p hybridiza- 
tion in the NiO lattice as compared to the Ni06 cluster. 
Anisimov at al. [H gave a similar argument to improve the 
agreement of their LDA+U calculation with the cluster 
calculations and in principle this should occur in any ap- 



proximation where a k independent self-energy is deter- 
mined in an 'impurity'-like calculation. More accurate 
knowledge about the total weight and width of the satel- 
lite probably will be needed to decide which is the more 
correct description. Finally it should be noted that the 
values of U and A determine mainly the energies of the 
upper Hubbard band and satellite relative to the valence 
band top - the band structure above —7eV is influenced 
hardly at all by these parameters. 

In summary: the variational cluster approximation due 
to Potthoff allows to combine the powerful cluster or CI 
method for transition metal compounds with the field- 
theoretical work of Luttinger and Ward to implement 
a variational scheme for the electronic self-energy and 
construct a band structure method for strongly corre- 
lated electron compounds. Both, a realistic band struc- 
ture and the full atomic multiplet interaction can be in- 
corporated into the Hamiltonian without problems, the 
system can be studied at arbitrarily low temperatures 
and the Green's function be obtained with arbitrary en- 
ergy resolution. The results are quite encouraging in that 
there is an essentially one-to-one correspondence between 
calculated Green's function and electron spectroscopies 
whereby a comparison of the fine structure of the broad 
band complex at the valence band top to experimental 
data of higher resolution would be desirable. The good 
agreement also suggests that the band structure of NiO is 
'Coulomb generated' in that the atomic multiplet struc- 
ture survives with minor modifications and gets broad- 
ened into weakly dispersive bands. All in all the VCA 
appears to be a promising tool for the study of transition 
metal compounds. 

The author would like to thank M. Potthoff for instruc- 
tive discussions. 
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